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Abstract 



In this letter we consider Andreev tunneling between a normal metal and a 
one dimensional Josephson junction array with finite-range Coulomb energy. The 
I — V characteristics strongly deviate from the classical linear Andreev current. We 
show that the non linear conductance possesses interesting scaling behavior when the 
chain undergoes a T = superconductor-insulator transition of Kosterlitz-Thouless- 
Berezinskii type. When the chain has quasi-long range order, the low lying excitation 
are gapless and the I — V curves are power-law (the linear relation is recovered 
when charging energy can be disregarded). When the chain is in the insulating 
phase the Andreev current is blocked at a threshold which is proportional to the 
inverse correlation length in the chain (much lower than the Coulomb gap) and 
which vanishes at the transition point. 



PACS numbers: 74.50, 74.60.Ge, 74.65+n 
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In the past few years there has been a renewed interest in the properties of meta- 
superconductor (NS) interfaces IQ. At low temperature and voltage the transport 
through a NS interface is due to Andreev tunneling 0]: a quasi-electron in the normal 
metal is reflected as a quasi-hole and a Cooper pair is added to the condensate on 
the superconducting side. Interference effects have been predicted for the subgap 
conductivity in NS tunnel junctions 0. 

Kim and Wen studied the effect of collective modes below the superconducting 
gap on the NS tunneling ||^; due to the excitation of the low-lying modes in the 
tunneling process, the I — V characteristics become non-linear. If the electrostatic 
energy is large enough, charging effects suppress transport at low temperature and 
voltage: this is the so called Coulomb blockade (see Ref.p] for reviews). Tunneling 
can take place only if electromagnetic modes with energy of the order of charging en- 
ergy can be excited whose effect on Andreev tunneling has been studied by Hekking 
et al. [0. 

In this letter we will study Andreev tunneling from a normal metal electrode to a 
one dimensional Josephson Junction Array (JJA). The system is depicted in Fig. 1. 
The transport in the JJA is due to Josephson tunneling. The Josephson energy 
(oc Ej) is minimized for a phase coherent superconducting state. If charging effects 
are relevant, the physics at low temperatures is dominated by its competition with 
the electrostatic energy {Eq = e^/2Co where Cq is the capacitance of the islands 
to the ground). The latter favours quantum fluctuations of the phase ip of the 
superconducting order parameter which can destroy superconductivity 0. In one 
dimensional JJA (which the case we will consider) this results in a T = phase 



transition p], ^ of the Kosterlitz-Thouless-Berezinskii type at ^Ej/8Eq = 2/tc. 

We will show that Andreev tunneling into a JJA is a powerful tool to perform 
spectroscopy of low lying modes of JJA. Moreover, it offers another way to study the 
critical properties of strongly interacting superconducting systems. Infact, quantum 
fluctuations of the phases in the JJA produce a nonlinear Andreev current which is 
very sensitive to the critical properties of the JJA. Here we point out some features 
of the phase diagram of the device which are reflected in the two-particle tunneling. 

Very recently it has been proposed that two electron tunneling can probe the 
spectrum of an 5 — 5* junction [|T^. The physics of this system is contained also in 
our setup in the limit of very small Ej. In this limit, the correlation length ^ of the 
JJA will not exceed few lattice constants and the transport properties are mainly 
determined by the junctions closest to the N-S boundary. 

We start from the Hamiltonian: 

- Ej^cos((/.,-¥.,+i)+t p^«7/;t(0)V^l(0) + /i.c.l . (1) 



It includes the Hamiltonians of the metal particle, that of the one dimensional 
JJA and the tunneling Hamiltonian which describes Andreev processes between the 
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normal metal and the first superconducting island respectively. In Eq. (|1|) ipl are 
the creation operators for the electrons of the normal metal particle and e(— iV) = 
— /i is the relative kinetic energy , measured from the Fermi energy. The integral 
is extended to the normal metal and a labels the spin index. Qi and cpi denote the 
charge and the phase of the i-th island of the JJA and are canonically conjugated, 
[Qi, ipj] = 2eSiji. We will consider a model for the capacitance matrix Cij where 
only the ground capacitance Cq and the nearest neighbor capacitance C are present. 
The range of the electrostatic interaction between Cooper pairs is A = ^C/Cq. In 
Eq.(|l|) tunneling occurs with an amplitude t at a single point of the junction |[TT[] . 



Finally we disregard quasiparticle tunneling in the JJA since we are dealing with 
energies far below the superconducting gap. 

The Andreev current across the NS boundary when a voltage V is applied is 
calculated from the time variation of the electron charge at the metal island 

Ij^g = -4e I t p ImDR{-2eV) 

/oo 
dte-^''^UQ{t) {[A^{t),A{0)]) (2) 
-oo 

where A(t) = e^'^'^^''^^ ip^{0,t)ipi{0,t). The retarded correlator in Eq. @ can be ob- 
tained from the analytic continuation of the corresponding time ordered Green 
function in imaginary frequencies. The latter is the Fourier transform of Dt{t) 
= {%A^iT)A{0)) = (r.^t(0,r)^}(0,0)) (r.V'i(0,r)^{(0,0)) ^(r) . All the key 
features are contained in the phase-phase correlation function of the first super- 
conducting island, g{T) = (7^e*'^°^^''e~*'^°''°'*). The critical properties of the JJA are 
reflected in the form of g{T) and can be studied by mapping the "quantum" one 
dimensional JJA at T = into a classical 1 + 1 dimensional (space and time) XY- 
model p] where Eq plays the role of the "temperature" . Vortices in the classical XY 
model correspond to T = phase slips in the JJA which occur at a certain junction 
of the chain and are due to quantum fluctuations modulated by Eq. 

Eq ^ Ej In this limit phase slips destroy the phase order and the Josephson 
chain is insulating. In the disordered phase the correlator g{T) decays exponentially 
in imaginary time. A careful evaluation of g{T) close to the transition point has 



been done in Ref. [0. At large times, taking into account the mapping we get 

g{r) = e-^'^'I^D/VK I r I)'/' r ^ oo (3) 



where Up = \J9>EjEq. In the limit C ^ Cq the correlation length is given by 



^ = exp[6/y2/7r — {Ej /SEq)^/'^] , where 6 ~ 1. It diverges at the superconductor- 
insulator transition. Inserting the analytic continuation of Eq.(§) into Eq.(0) we get 
the tunneling current at small voltages of the array in its insulating phase 

I^s = sign(y) (2e | V \ /cUp - C') [2e | V \ /up - C'f' . (4) 
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Here Ga = Svre^ | t ^ A^^(O) is the Andreev conductance for a bulk superconductor 
(A^(0) is the density of states of the metal particle). The typical I — V characteristics 
is shown in Fig. (2) , curve (a). The exponential decay of correlation in the JJA cor- 
responds to the Coulomb blockade of the Andreev tunneling. Interestingly enough, 
the threshold is not the bare charging energy but it is proportional to the inverse 
correlation length. Hence it vanishes at the transition with a behavior related to 
the critical state of the chain, 

As long as the correlation length ^ is larger than A, the I — V characteristics are 
mainly determined by the self-capacitance and eq.(^) is valid also for finite C. At 
lower values of the ratio Ej/Eq further away from the transition, the correlation 
length ^ decreases and eq.(|) does not hold any longer because C cannot be disre- 
garded. 



At the critical value yEj/SE^ = 2/tt the JJA undergoes a transition to the 
superconducting phase where the phase-slips of opposite sign are tightly bound in 
pairs and the voltage drop across the array is zero. 

Ej ^ Eq The Josephson Hamiltonian is well approximated by its phase- 
wave Gaussian limit in which quantum phase slips are neglected. The low energy 
excitations are modes with dispersion 



2 ^ 2Ej(l-cos q) 



The resulting phase correlator in imaginary time is 



where f = KUpT, k = (1/4 + C/Coy^/^ and r/ = 1 + {2%)'^ ^SEq/Ej. The system 
possesses quasi-long range order and the correlations decay with a power-law for 
r ^ oo, gi^r) ~ | UpT \ ''^^ . As a result the I — V characteristics at low voltage have 
a power-law behavior 

^ e-^(^-i) Up f2e\V\y 



r(l +r]) 2e \ u 



p 



where 7 is the Euler number. In the classical limit, 77 — ^ 1 , the ground state of 
the chain has all the phases aligned and no collective modes are excited. Then the 
chain behaves as a bulk superconductor and the Andreev I — V characteristic is 
linear. When Eq increases phase-wave collective modes become active. In particular 
they are excited by the Andreev tunneling events themselves and their back-action 
determines the anomalous nonlinear characteristics of Eq. (|^. 
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The phase-wave Gaussian hmit is the fixed point Hamiltonian also when quantum 
phase shps are included but the parameters are renormalized. In Eq. (^) rj is replaced 
by the renormalized value which deviates from the phase- wave value close to the S-I 



transition (see Ref.||T3| for more details). 

Comparing eq.(^) with eq.(|^ one sees that close to the transition point the 
I — V curves start as a power law, but with different exponents when approaching 
the transition from the two sides (see the curves (b) and (c) in Fig. (2)). Approaching 
the transition from above r/^ = 5/4. The jump is 

Vcr-Vtr = 1/4 

The jump the I—V curves is related to the jump in the superfluid density at the tran- 
sition point. A similar behavior, though with a different value, has been predicted 
and observed long time ago for the conductance in the classical 2D superconducting 
films 0. 

Now we analize in more detail the Andreev I — V characteristics in the ordered 



phase, ^JEj/^iEq > 2/7r. Eq. (|^ is vahd at voltages much smaller than uip/e. For 
voltages much larger than Up/n the short time limit of Eq.(^) is needed, which is 

^?(r)=exp{-^(7rf-f2 + iL^3)| 
Analytic continuation leads to 

TV , ^ , KiOr, \ rr , KCUp 



lNs = G^iV--iv-l)^ (9) 
V 2 e / e 

with corrections which vanish exponentially with the applied voltage. At very large 
voltages, where eq. applies, we recover the linear behavior of the I — V character- 
istics, apart from an offset which is related to the average energy cost for exciting 
the collective phase-wave modes of the chain. The I — V curves are shown in the 
inset of Fig. (2). We calculated the behavior at intermediate voltages directly in real 
times from Eq.(^. The excitations of the chain act as bosonic modes of an electro- 
magnetic environment |jT5[ coupled to tunneling pair and their energy spectrum is 



given by Eq (^ . Integrating the Bose operators of the environment we get the real 
time correlator 

P{t) = {e'^o{t)^-iMo)^ = exp j-l" ^ C(g) [1 - cosiiu.t) + tsm{ugt)]\ (10) 



where ({q) = ye^C~^(g)/i?jsin^(g/2). The current takes the form 

'•2ev duj 



l-NS 



nzev 

-4e\t\^N\0) — {2eV - u) P{u) . (11) 

Jo 271 



The power law behavior of Eq. (0) for small voltages turns very rapidly in the asymp- 
totic form, Eq. (^, at ~ nUp/e. P{uj) can be interpreted as the probability that 
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a tunneling event excites a mode of the JJA. This results resembles that obtained in 
connection with the study of the effect of the electromagnetic environment on single 



electron tunneling [|T5|. In this approximation the one-dimensional JJA acts like a 

dance 

e{Kujp-uj) . (12) 



line of linear elements with the equivalent impedance 

Z(uj) V - 1 



In the phase- wave approximation the function P{uj) presents an elastic peak, P(u;) ~ 
uj^~'^ for C(j — i> 0, and a divergence at u = nup, P{uj) ~ \Kujp — uj\'^~^^'^. This latter 
feature is associated with excitations which are almost localized in some junction 
of the JJA and whose density of states diverges for the infinite chain. It produces 
some structure near V = nup/e both in d^I/d^V ~ P{eV) and in the differential 
conductance dl/dV which is shown in the inset of Fig. (3). Notice that the low 
voltage behaviour is dominated by the capacitance to the ground Co since C enters 
simply in the scaling factor k (see Eqs.(|9|, p^) ). 

The modifications we expect in going beyond the phase-wave approximations in 
the ordered state can be qualitatively discussed in terms of the spectrum of the exci- 
tations we neglected. Creation and breathing of bound pairs of phase slips produce 
low energy Gaussian fluctuations which, as already discussed, only renormalize t] 
and modify the exponent in Eq. (|^) accordingly. They should not affect the struc- 
ture near V = nup/e which is due to localized excitations. The latter, however, 
can be modified due to the interaction of these modes with the phase slips or be- 
tween themselves. Moreover going to higher voltages phase slips may be induced. 
Their typical energy scale, away from the transition, is larger than nup so in this 
regime source induced phase slips are believed not to be important. They should 
start to affect the behavior at ~ KUp/ e only close to the transition, so we expect 
a detectable region in the phase diagram where the behavior of Eq. (^) can be 
observed. At even higher voltage all the above excitations can be active and this 
should increase the offset in the I — V curve, Eq. (^). 

In conclusion we have shown that Andreev tunneling may be a powerful method 
to investigate the phase diagram and the excitation spectrum (see also [|l^) of JJA. 
The results derived in this letter are summarized in Figs. (2, 3). When the JJA is 
in the disordered phase the Andreev tunneling is blocked up to a threshold voltage 
proportional to the correlation length ^ (Fig. (2)). At the superconductor- insulator 
transition the I — V curves are power law like with a jump in the exponent at 
the transtion point (Fig. (2)). When the JJA is in the superconducting phase, 
the characteristics start as a power law and approach rapidly the shifted ohmic 
behaviour. A richer structure generated by localized excitations of the JJA is seen 
in the derivatives of the I—V curves (see Fig. (3)), near V = nujp/e. We have studied 
explicitely a one dimensional JJA but other types of arrays can be treated along 
the same lines. In two dimensions frustration effects due to an applied magnetic 
field play an important role. The interference patterns manifest in the differential 
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conductance in a rather peculiar way. Work in this direction is in progress. 
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FIGURE CAPTIONS 



Figure l:The system considered in the present work: a normal metal island con- 
nected by a tunnel junction to a very long chain of Josephson junctions. 

Figure 2: (a) The I -V curve in the limit Eq '> Ej: I = ioi V < Vth- At the 
transition point the current starts as a power law, however there is a jump in the 
exponent of the characteristics (curves (b) and (c)). For £'j — > oo the classical linear 
behavior is recovered (d). In the inset some curves in the Ej ^ Eq regime. 

Figure 3:The function P{oj), which also gives the second derivative of the I — V 
curves, in the limit Ej > Eq. Eq/Ej = 0.05, 0.02, 0.005, 0.0005, from top to bot- 
tom. In the inset the differential conductance, in the same limit, is suppressed for 
increasing Eq/ Ej. 
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